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Abstract 

The Navier-Stokes equations are considered by the use of the method 
of Lagrangians with covariant derivatives (MLCD) over spaces with affinc 
connections and metrics. It is shown that the Euler-Lagrange equa- 
tions appear as sufficient conditions for the existence of solutions of the 
Navier-Stokes equations over (pseudo) Euclidean and (pseudo) Rieman- 
nian spaces without torsion. By means of the corresponding (n — 1) + 1 
projective formalism the Navier-Stokes equations for radial and tangential 
accelerations are found. 

1 Introduction 

By the use of the method of Lagrangians with covariant derivatives (MLCD) pQ, 
the different energy-momentum tensors and the covariant Noether's identities 
for a field theory as well as for a theory of continuous media can be found. 
On the basis of the (n — 1) + 1 projective formalism and by the use of the 
notion of covariant divergency of a tensor of second rank the corresponding 
covariant divergencies of the energy-momentum tensors could be found. They 
lead to Navier-Stokes' identity and to the corresponding generalized Navier- 
Stokes' equations. 

The general scheme for obtaining the Navier-Stokes equations for radial and 
tangential acceleration could be given in the form 
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The structure of a Lagrangian theory of tensor fields over a differentiable 
manifold M (dim M = n) could be represented in the form 
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Let the following structure 

(M, V, g, T, P) 

be given, where 

(i) M is a differantiable manifold with dim M = n, 

(ii) V = V A b ■ eA®e B G <8> fc j(M) are tensor fields with contravariant rank 
k and covariant rank I over M, A and B are collective indices, 

(iii) g € < 8>Sym2(M) is a covariant symmetric metric tensor field over M, 

(iv) r is a contravariant affine connection, P is a covariant affine connection 
related to the covariant differential operator along a basis vector field <9j or 
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in a co-ordinate or non-co-ordinate basis respectively 



V dt V = V A B , t -d A ®dx B 

V A B;i = (V A B ,i +T A r V C B + P° ■ V A D 



V A 



_ dV A j 
B ' 1 dx 1 



A Lagrangian density L can be considered in two different ways as a tensor 
density of rank with the weight q = 1/2, depending on tensor field's compo- 
nents and their first and second covariant derivatives 

(i) As a tensor density L of type 1, depending on tensor field's components, 
their first (and second) partial derivatives, (and the components of contravariant 
and covariant affine connections), i.e. 



L — V '~dg ■ L(gij, gij,k,9ij,k,h V b, V B,i, V B,i,j, Ij fc , P l jkl 
where L is a Lagrangian invariant, 

d g = det(gij) < , g — gij ■ dx % .dx 3 , 
dx l .dx J = - ■ (dx 1 <g> dx j + dx 3 ® dx 1 ) , 



V A 



dV A 

B,i,j 



dx J dx 1 



The method using a Lagrangian density of type 1 is called Method of La- 
grangians with partial derivatives (MLPD). 

(ii) As a tensor density L of type 2, depending on tensor field's components 
and their first (and second) covariant derivatives, i.e. 



L — \J—d g ■ L(gij, gij;k, gij;k;U V B, V B;i, V B;i;j) ■ 

By the use of the variation operator 6, commuting with the covariant differ- 
ential operator 

6oVs = V ? o(5 + V 5c , eeT(M) , T{M) = U xeM T x (M) , 

we could find the Euler-Lagrange equations. 

By the use of the Lie variation operator (identical with the Lie differential 
operator) we could find the corresponding energy-momentum tensors. 

The method using a Lagrangian density of type 2 is called Method of La- 
grangians with covariant derivatives (MLCD). 

1.1 Euler-Lagrange's equations 

The Euler-Lagrange equations follow from the variation of the Lagrangian den- 
sity of type 2 in the form |H] 
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(i) for the tensor fields V 



P A B = 



SV A B 

(ii) for the metric tensor field g 
ogki 2 

Special cases: (Pseudo) Euclidean and (pseudo) Riemannian spaces without 
torsion. 

S « L - n 5gL i 1 t o M -n 

5V* B -° ' 5g kl + 2 ~ ■ 

1.2 Energy-momentum tensors 

By the use of the Lie variation operator the energy-momentum tensors follow 

m 

(i) Generalized canonical energy- momentum tensor 9 = 9i 3 ■ dj <g> dx l , 

(ii) Symmetric energy- momentum tensor of Belinfante S T = Ti 3 ■ dj ® dx l , 

(iii) Variational energy-momentum tensor of Euler-Lagrange Q — Q i 3 ■ dj <8> 

The energy-momentum tensors obey the covariant Noether identities 

T, ■ T), I) , 

F + 69 = , 

(first covariant Noether's identity) 

9% ^ 3 = 3 j 

9- S T = Q . 

(second covariant Noether's identity) 

Now we can draw a rough scheme of the main structure of a Lagrangian 
theory: 

n 



r < — L — > 

i S T 6L/6V A B 

9 | Q 

l — ► 9- S T = Q <—- 

I i 

L — > F + S9 = < — j 
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2 Invariant projections of the energy- momentum 
tensors 



By the use of the (n — 1) + 1 projective formalism we can find the invariant 
projections of the energy-momentum tensors corresponding to a Lagrangian 
field theory or to a theory of continuous media. The idea of the projective 
formalism is the representation of the dynamic characteristics of a Lagrangian 
system by means of their projections along the world line of an observer and to 
local neighborhoods orthogonal to this world line. The tangent vector to the 
world line of the observer and its local neighborhoods determine the notion of 
frame of reference [5] Fr(u, r, where u is the tangent vector of the world 
line, r is the parameter of the world line, interpreted as the proper time of the 
observer £j_ is a contravariant vector field, orthogonal to u. The variation of £j_ 
determines the relative velocity and the relative acceleration between the points 
at the world line and the points in the neighborhoods lying in the sub space 
orthogonal to the vector u. 

Let the contravariant vector field u s T(M), g(u,u) := e / 0, and its 
corresponding projective metrics h u and h u 

h u = g • g{u) ® g(u) , h u = ~g ■ u® u , 

e e 

g = g l] -d i .d j , di-dj = - • (<% <8> dj + dj ® <%) , 

be given. Then the following proposition can be proved: 

Proposition 1. Every energy-momentum tensor G ~ (6, S T, Q) could be 
represented in the form [3] 

G = (pa H L ■ k) ■ it® g{u) - L-Kr + u® g( k n) + k s <E> g(u) + ( k S)g , 

e 

where 

G— G— kn G~q 

TT = TT , S = S, D= D , 

Pg is the rest mass density, k = (1/e) • [g(u)](Kr)u, L is the pressure of 
the system, Kr — g l . ■ di <E> dx 3 is the Kronecker tensor, k n is the conductive 
momentum density, k s is the conductive energy flux density, k S is the stress 
tensor 0]. 

3 Covariant divergency of the energy- momentum 
tensors and the rest mass density 

The covariant divergency SG of the energy-momentum tensor SG (G ~ 9, S T, 
Q) can be represented by the use of the projective metrics h u , h u of the con- 
travariant vector field u and the rest mass density for the corresponding energy- 
momentum tensor pc- 
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8G and g(SG) can be found in the forms [3] 

SG= { PG + \-L-k)-g{a)+ 
+Hp G + i ■ L ■ k) + {p G + ±-L-k)-Su + 5 G s}- g(u)- 
- KrL - L ■ SKr + 5u ■ g( G W) + g(V u G W) + g(Va-gu)+ (1) 
+ {pa + \ ■ L ■ k) L (Vu5)(«) + (V uff )( G 7f) + {Vg- s9 ){u) + 
+ 5(( G S)g) , a = V„u , 

g(SG) = (p G + i • i • k) ■ a+ 
+ [u(p G + \-L-k)+ ( PG + ±-L-k)-Su + 6 G s]-u- 
- g( KrL) - L ■ g{8Kr) + Su ■ % + V„ G n + V Gl rU+ (2) 
+ (p G + \ ■ L ■ k) ■ g{V u g){u) +j(V uff )( G ^) + g{V G ^g)(u) + 
+ g(5(( G S)g)) ■ 
In a co-ordinate basis SG and ~g(SG) will have the forms 

G l 3 d = (p G + \ ■ L ■ k) ■ a t + 
+ [(pa + \-L- k) d -ui + {p G + \-L-k)- v? d + G s J -j] ■ Ui - 
-L A - L ■ glj + ;j ■ G lf l +g t j ■ ( G 7fJ , k _- u k + vP jh ■ G s k ) + (3) 
+ 9ij;k • [(PG + \-L-k)-u' -u k + G 7f J • u k + v? ■ G s k ]+ 
+ (9ik- G S jk h, 

f k Gt i- . . = ( PG + ±.L-k)-a*+ 
+ [{pa + \ ■ L ■ k) d ■ u> + {p G + \ ■ L ■ k) ■ ui .j + G s> , 3 ] ■ u *- 
—Lj ■ g*i -L-g lk ■ g{ ;j + v? ;j ■ G 7f l + g tt z ;j ■ u j + u l .j ■ G s j + (4) 
+ 9 a ■ giy.k ■ [{pa + \ ■ L ■ k) ■ u~ ■ u k + G # • u k + vJ ■ G s fc ] + 

+ g fl -(gr k - G s>% ■ 

4 Navier-Stokes' identities and Navier-Stokes' 
equations 

If we consider the projections of the first Noether identity along a non-null 
(non-isotropic) vector field u and its corresponding contravariant and covariant 
projective metrics h u and h u we will find the first and second Navier-Stokes 
identities. 

From the Noether identities in the form 

~g(F) +~g(58) = 0, (first covariant Noether's identity) , 

(6)g — ( s T)~g = (second covariant Noether's identity) , 

we can find the projections of the first Noether identity along a contravariant 
non-null vector field u — u l ■ di and orthogonal to u. 
Since 

g(g(F) , u) = g lk ■ g Yl ■ F x ■ u l = g\ ■ F, • u 1 = F % ■ u 1 = F(u) , (5) 
g(g(59),u) = (66)(u) , F = F k -dx k , (6) 
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we obtain the first Navier-Stokes identity in the form 

F(u) + (66)(u) = . (7) 

By the use of the relation 

g[h u (g)(F)]=g(h u [g(F)}) = h u (F) , g(h u )g = h u , (8) 
g[h u (g)(S9)} = g(h u [g(69)}) = h«(50) , (9) 

the first Noether identity could be written in the forms 

h u [g{F)]+h u [g{56)] =0 , (10) 
h u {F) + h u {56) =0 . (11) 



The last two forms of the first Noether identity represent the second Navier- 
Stokes identity. 

If the projection h u (F), orthogonal to u, of the volume force F is equal to 
zero, we obtain the generalized Navier-Stokes equation in the form 

h u {56) = , (12) 

or in the form 

h u [g{56)}=0 . (13) 

Let us now find the explicit form of the first and second Navier-Stokes iden- 
tities and the explicit form of the generalized Navier-Stokes equation. For this 
purpose we can use the explicit form of the covariant divergency 89 of the gen- 
eralized canonical energy-momentum tensor 9. 

(a) The first Navier-Stokes identity follows in the form 

F(u) + {p e + - ■ L ■ k) ■ g(a,u)+ 
e 

+ e ■ [u(p e + -■ L ■ k) + (p e + - ■ L ■ k) ■ 5u + 5 e s}- 

e e 

- (KrL)(u) - L ■ (SKr)(u) + g(V u 9 tF, u) + g(Ve^u, u)+ 

+ (p e + --L-k)- (V u g)(u,u) + (V u g)( 8 Tf,u) + (Vo^g)(u,u)+ 
e 

+ [S(( e S)g)}(u)^0 . (14) 

Since 

g(u, a) = ±l u ■ ^ - ^ • (S7 u g)(u, u) = 
= \-[^{±ii)-{V u g){u,u)] , 

the first Navier-Stokes identity could be interpreted as a definition for the change 
of l\ along the world line of the observer. The length of the non-isotropic con- 
travariant vector u is interpreted as the velocity of a signal emitted or received 
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by the observer [S]. On this basis, the first Navier-Stokes identity is related to 
the change of the velocity of signals emitted or received by an observer moving 
in a continuous media or in a fluid. 

(b) The second Navier-Stokes identity can be found in the form 

h u [g(F)} + h u [g(S9)} = 



= (pe + - ■ L ■ k) ■ h u (a)- 

e 

- h u [g(KrL)] - L ■ h u [g(SKr)} + Su ■ h u ( 6 Tf)+ 

+ Mv u e 7f) + K(y •»- s u)+ 

+ (p e + - ■ L ■ k) ■ h u [g{V u g)(u)} + h u [g(\7 u g)( 9 n)} + 

e 

+ h u [g(V^g)(u)] + h u [g(K( e S)g))} + K[g{F)] = . (15) 

(c) The generalized Navier-Stokes equation h u \g(50)] = follows from the 
second Navier-Stokes identity under the condition h u [g(F)] = or under the 
condition F = 

(pe + - ■ L ■ k) ■ h u (a)- 
e 

- h u [g{KrL)\ - L ■ h u [g{5Kr)\ + Su ■ h u ( e W)+ 
+ h u (V ' u e Tf) + h u {V , e- s u) + 

+ (p e + --L-k)- h u [g{V u g)(u)} + h u {g(\7 u g)( 9 n)} + 
e 

+ h u [g(V^g)(u)} + h u [g(5(( 6 S)g))} 



= , (16) 

hu{a) =g(a) g{u,a) ■ g(u) . (17) 

e 

The second Navier-Stokes identity could be considered as a definition for the 
density of the inner force. If the density of the inner force is equal to zero, i.e. 
if F = g{F) = 0, then the covariant divergency, 59 — g{89) of the generalized 
canonical energy-momentum tensor is also equal to zero, i.e. 69 = g(59) = 0. 
Then the orthogonal to the contravariant vector field u projection of the second 
Navier-Stokes identity lead to the equations 

g[h u (F)}=0 &g[h u (59)}=0 . 

The last equation is the Navier-Stokes equation in spaces with affine connec- 
tions and metrics. Now, we can prove the following proposition: 

Proposition 2. The necessary and sufficient condition for the existence of 
the Navier-Stokes equation in a space with affine connections and metrics is the 
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condition for the vanishing of the density of the inner force in a dynamic system 
described by the use of a Lagrangian invariant L, interpreted as the pressure p 
of the system, i.e. the necessary and sufficient condition for 

g[h u (66)] = 

is the condition 

g[h u (F)] = 

The proof follows directly from the projective second Navier-Stokes identity 
g[h u (F)}+g[h u (Se)}=0. 

Special case: (L n , g)-spaces: S ~ C, f l j — g % -, g{u,u) = e = const. ^ 0, 
k = l. 

SKr = , (18) 

(a) First Navier-Stokes' identity 

F(u) + (p g + --L)-g{a,u)+ 

e 

+ e ■ [u(p e + - ■L) + (p g + --L)-5u + 5 e s}- 
e e 

- (KrL)(u) + g(V u °tt, u) + g{Ve^u, u) + 

+ { Pe + --L)- {V u g)(u, u) + {V u g){ d n, u) + (V»,ff)(u, u)+ 
e 

+ [S(( e S)g)}(u)^0 . (19) 

(b) Second Navier-Stokes' identity 



(po + - ■ L) ■ h u (a)- 
e 

- h u [g(KrL)] + Su ■ h u { e if)+ 
+ h u (V u e W) + h u {Ve^u)+ 

+ {pe + \-L)- h u [g{V u g){u)] + h u [g{V u g)( e n)} + 
+ h u [g(Vs 1 g)(u)} + h u [g(S(( e S)g))} 



+h u [g(F)}=0. (20) 
(c) Generalized Navier-Stokes' equation h u \g(59)\ = 
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[pe + - ■ L) ■ K(a)- 

- h u [g(KrL)} + Su ■ h u ( e Jr)+ 
+ h u (X7 u 9 Tf) + h u (Ve- s u)+ 

+ {p e + j ■ L) ■ h u [g{\7 u g)(u)] + h u [g{V u g)( e W)} + 
+ h u [g(Vs^g)(u)} + h u [g(S(( 9 S)g))} 



= . (21) 

Special case: V^-spaccs: S = C, f % j = gj, \7^g = for E T(M), g(u,u) = 
e = const, 0, k = 1, g(a, u) = 0. 
(a) First Navier-Stokes' identity 



F{u) + 

+ e ■ [u(p e + - ■ L) + (p e + - ■ L) ■ 5u + 5 e s] - (KrL){u)+ 

+ [S(( e S)g)}(u)^0 . (22) 

(b) Second Navier-Stokes' identity 

{pe + \-L)- h u (a) - h u [g(KrL)} +8u-h u ( e W) + h u [g(S(( 8 S)g))} + K [g(F)] = . 

(23) 

(c) Generalized Navier-Stokes' equation h u [g(59)\ = 

(pe + \-L)- h u {a) - h u [g(KrL)] + 5u ■ h u ( 6 W) + h u [g(S(( e S)g))} = . (24) 

If we express the stress (tension) tensor ( e S)g by the use of the shear stress 
tensor fc s D, rotation (vortex) stress tensor ', and the expansion stress in- 
variant kU then the covariant divergency of the corresponding tensors could be 
found and at the end we will have the explicit form of the Navier-Stokes identi- 
ties and the generalized Navier-Stokes' equation including all necessary tensors 
for further applications. The way of obtaining the Navier-Stokes equations could 
be given in the following rough scheme 
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5 Invariant projections of Navier-Stokes' equa- 
tions 

5.1 Navier-Stokes' equations and Euler-Lagrange's equa- 
tions 

Let us now consider the second Navier-Stokes identity in the form |3] 
g[h u [g(F)]]+g[h u [g(8e)}]=0 

or in the form 

F± + 60 x = , 

F x =g[h u [g(F)]] , 56 x = g[h u [g(S0)} , 
g(u,F ± ) = , g{u,S0±)=O . 

The explicit form of the density F of the inner force could be given as [3] 

F = Fi- dx i , 

Wi - Wi{T m \ g jk -,i) , 
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where TJy 3 are the components of the torsion tensor (in a co-ordinate basis 

r« J ' = r/ fc -rj,). 

For (pseudo) Euclidean and (pseudo) Riemannian spaces without torsion 
(T^ = 0) the quantity W is equal to zero (Wi = 0) and the density of the inner 
force F has the form 



5L 



rA 



If the Euler-Lagrange equations are fulfilled in (pseudo) Euclidean and (pseudo) 
Riemannian spaces without torsion, i.e. if 

= , 



SV A L 
then F = and 

F i = ^-.V- B , = , 

and the following propositions could be proved: 

Proposition 3. Sufficient conditions for the existence of the Navier-Stokes 
equation in (pseudo) Euclidean and (pseudo) Riemannian spaces without torsion 
are the Euler-Lagrange equations. 

Proposition 4- Every contravariant vector field u G T(M) in (pseudo) Eu- 
clidean and (pseudo) Riemannian spaces without torsion is a solution of the 
Navier-Stokes equation if the Euler-Lagrange equations are fulfilled for the dy- 
namic system, described by a given Lagrangian invariant L = p interpreted as 
the pressure of the system. 

Corollary. If L = p = p(u\ u l -j, u l ;J;fe , g tj , g lj]k , Qij-k-,1, V A B , V A B ,i, 
y A B-,i;j) is a Lagrangian density fulfilling the Euler-Lagrange equations for 
u l and V A b in (pseudo) Euclidean and (pseudo) Riemannian spaces without 
torsion, then the contravariant non-isotropic vector field u is also a solution of 
the Navier-Stokes equation. 



5.2 Representation of F± and 59 ± 

Now, we can use the corresponding to a vector field g(u, = (orthogonal 
to the vector field u) projective metrics ti£ ± and h^ 1 - 

kx=9- ( , t , ■ g($±) ® 9 (gx) > 

h u =9- u * s ■ U <8 U ■ 
The vector field F± could be written in the form [J, [S] 
F ± = 4^44 • ^ + 9[hu (F±)] = T9(Fx,n ± ) ■ n± + g[h u {F ± )\ = 
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F± z is the radial inner force density and F± c is the tangential (Coriolis) 
inner force density 

F± z = T9{F±,n±) ■ n± , F ±c = g[h u (F±)] , 
g(F ±z ,u) = , g(F ±c ,£ ± ) = , g(F ±c ,u) = . 

The Navier-Stokes equation could now be written in the form 

56a_ = T9(69±,n±) ■ n± + g[h u (59 ± )] = , 

or in the forms 

50j_ z := Tg(S0j_, n±_) ■ n± = , 

Navier-Stokes' equation for radial accelerations , 

50 ±c :=g[h u (50±)]=O , 

Navier-Stokes' equation for tangential accelerations. 

5.3 Radial projections of Navier-Stokes' equation. Navier- 
Stokes' equation for radial accelerations 

If we use the explicit form of the Navier-Stokes equation 
(pe + ^ • L ■ k) ■ a±- 

- \g(KrL)] ± - L ■ [g(6Kr)} ± + 5u ■ e ir ± + 
+ (V u V)± + (V»*u)x+ 

+ (p e + - • L ■ k) ■ [g(V u g)(u)} ± + [g(V u g)( e n)} ± + 
e 

+ [9(V°s9)(u)} ± + [g(S(( e S)9))}± = , (25) 
L=p, (26) 

and apply the projection of the Navier-Stokes equation along and orthogonal to 
the vector field by the use of the representation of the acceleration a± in 
the form 

a± = g(a±,n±) ■ n ± +g[h u (a ± )} = a z + a c , 
a z = g(a±,n±)-n± , a c =g[h^ ± (a±)] , 

where a z = g(a±,n±) ■ n± = z fl az • n± is the radial (centrifugal, centripetal) 
acceleration and a c = g[h^ ± (a±)} = Tla c ■ m±, g(n±,m±) = 0, is the tangen- 
tial (Coriolis) acceleration, we could fine the explicit form of the Navier-Stokes 
equation for radial (centrifugal, centripetal) accelerations in the form 
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(pe + - ■ L ■ k) ■ a z - 

e 

~ [g(KrL)] ±z - L ■ [g(SKr)] ±z + 5u ■ e n ±z + 
+ (V„ 9 7f) ±2 + (VosU)± z + 

+ (p e + --L-k)- [g(V u g)(u)] ±z + [g(V u g)( e w)] ±z + 

e 

+ &(V^5)(«)]x, + [g(S(CS)g))h z = , (27) 
L = P , (28) 

Special case: Perfect fluids: % = 0, e s = 0, e S = 0, L = p. 



(p g + --L-k)-a z - [g(KrL)} ±z - L ■ [g(6Kr)} ±z + 

e 

+ ( Pe + - ■ L ■ k) ■ [g(V u g){u)]xz 
e 

= . (29) 

Special case: Perfect fluids in (pseudo) Euclidean and (pseudo) Riemannian 
spaces without torsion: e W = 0, e s = 0, e S = 0, L = p, \7 u g = 0, SKr = 0. 

(pe + - ■ p) ■ a z = [g(Krp)]± z , 
e 

[g(Krp)]± z . 



5.4 Tangential projections of Navier-Stokes' equation. Navier- 
Stokes' equation for tangential accelerations 

For tangential (Coriolis') accelerations the Navier-Stokes equation takes the 
form 0] 

(pe + - ■ L ■ k) ■ a c - 

e 

- [g(KrL)]tc - L ■ [g(5Kr)] ±c + Su ■ e n^ c + 
+ (V„ 9 7f)±c + (V^u)±c+ 

+ {pe + - ■ L ■ k) ■ [g(V u g)(u)]xc + [ff(V u5 )( e 7f)] ±c + 

+ [9(Ves9)(u)]xc + [g(S(( 6 S)g))hc = , (30) 
L=p. (31) 
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Special case: Perfect fluids: °n = 0, s = 0, S = 0, L = p. 



[pe + ^ • L ■ k) ■ a c - 

- [g(KrL)] ±c - L ■ [g(5Kr)] ±c + (p e + - ■ L ■ k) ■ [g(V u g)(u)} ±c = . (32) 

e 

Special case: Perfect fluids in (pseudo) Euclidean and (pseudo) Riemannian 
spaces without torsion: 7f = 0, s = 0, e S = 0, L = p, V u g = 0, SKr = 0. 



(pe + - ■ p) ■ a c = 

e 

= [g(Krp)] ±c , (33) 
1 



[pe + --p) 



■[g(Krp)] ±c . (34) 



6 Conclusions 



The representations of the Navier-Stokes equation in its forms for radial (cen- 
trifugal, centripetal) and tangential (Coriolis') accelerations could be used for 
description of different motions of fluids and continuous media in continuous 
media mechanics, in hydrodynamics and in astrophysics. The method of La- 
grangians with covariant derivatives (MLCD) appears to be a fruitful tool for 
working out the theory of continuous media mechanics and the theory of flu- 
ids in spaces with affine connections and metrics, considered as mathematical 
models of space-time. 
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